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ABSTRACT 

The  objective  of  the  thesis  is  to  minimize  the  total 
time  required  for  a  driver  to  pass  through  the  Lincoln 
Tunnel,  travel  to  his  destination  in  the  Central  Business 
District  of  New  York  City,  wait,  be  serviced  and  then  exit 
his  destination.   This  is  accomplished  by  requiring  vehicles 
to  delay  in  New  Jersey  before  entering  the  tunnel. 

A  general  model  for  total  time  in  the  system  as  a 
function  of  the  waiting  time  in  New  Jersey  is  presented.   A 
sampling  technique  is  developed  which  tests  the  adequacy  of 
that  portion  of  model  which  predicts  the  conditions  within 
New  York  City. 

The  control  techniques  are  extended  to  apply  to 
ships  entering  harbor  through  the  constrains  of  a  channel. 
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CHAPTER  1 
DESCRIPTION  OF  THE  SYSTEM 

I. 1.   Introduction 

The  central  business  district  (C.B.D.)  of  New  York 
City,  south  from  42nd  Street,  may  be  viewed  as  a  large  and 
complex  commerical  and  industrial  center.  Over  the  years, 
activities  of  a  similar  nature  have  tended  to  locate 
together  and  thus  somewhat  separate,  although  not  entirely 
independent,  business  districts  have  formed. 

There  are  numerous  entrances  into  and  exits  out  of 
the  C.B.D.   From  the  point  of  view  of  the  vehicles  travel- 
ing into  New  York  City,  the  various  business  districts  may 
be  thought  of  as  traffic  receptors  or  traffic  sinks.   Hence- 
forth in  this  thesis,  these  will  be  referred  to  as  traffic 
sinks  or  simply  as  sinks.   At  the  present  time  there  is  no 
attempt  made  to  limit  the  number  of  vehicles  that  enter  the 
C.B.D.  and,  as  a  result,  large  traffic  congestions  are 
common.   Consequently,  a  driver  entering  New  York  City  can 
expect  long  delays  before  concluding  his  business. 

The  attempt  of  this  thesis  is  to  reduce  these  delays 
by  requiring  vehicles  to  wait  before  entering  the  C.B.D. 
This  will  allow  the  particular  traffic  sink  to  clear  out 
enough  so  that  the  total  time  a  driver  spends  in  New  York 
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City  will  be  reduced.   One  entrance,  the  Lincoln  Tunnel, 
will  be  studied  in  detail. 

1.2.  Business  Districts  in  New  York  City 

Although  no  rigorous  attempt  has  been  made  to 
delineate  the  New  York  City  traffic  sinks,  the  Port  of  New 
York  Authority  has  supplied  the  author  with  sufficient  infor- 
mation to  describe  roughly  the  boundaries  as  shown  in  Fig.  1. 

The  author  attempted  to  determine  exactly  the  limits 
of  one  of  these  sinks.   It  was  found  that  while  difficult 
and  time  consuming,  such  an  undertaking  was  not  impossible 
and  could  be  effected  if  sufficient  personnel  and  time  were 
available.   Fig.  1  is  time  dependent;  for  example,  districts 
12  and  13  are  in  the  process  of  relocating  further  north. 
Other  districts  can  be  expected  to  expand,  contract  or 
relocate  as  time  goes  on;  this  indicates  continuous  updat- 
ing. 

1.3.  The  Lincoln  Tunnel 

Some  of  the  entrances  to  the  New  York  City  C.B.D. 
are  also  shown  in  Fig.  1.   The  Lincoln  Tunnel  is  a  three 
tube,  six  lane  facility  operated  by  the  Port  of  New  York 
Authority  and  represents  a  cost  of  approximately  $18,000,000 
per  lane  per  mile.^  '   £ach  tube  has  2  lanes  and  is  about 
7,650  feet  long.   The  south  tube  accommodates  east  bound 
traffic  exclusively,  the  north  tube  handles  only  west  bound 
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FIGURE  1 
The  Traffic  Sinks  of  New  York  City 
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traffic,  and  the  center  tube  is  switchable.  During  the 
morning  rush  hour  traffic,  the  center  tube  is  used  to  trans- 
port vehicles  into  New  York  and  during  the  evening  rush  hour 
it  is  changed  into  a  west  bound  roadway.  At  other  times  the 
central  tube  handles  traffic  flows  in  both  directions. 

Although  rather  complete  information  on  the  origins 
and  destinations  of  vehicles  which  use  the  Lincoln  Tunnel 
has  been  collected^")  and  data  on  the  behavior  of  vehicles 
within  the  tunnel  has  been  obtained, (26)  no  8UCCQSS  has  been 
realized  in  attempting  to  optimize  the  use  of  the  center 
tube/11) 

Since  the  adequacy  of  the  control  techniques,  to 
be  presented  in  later  chapters,  depends  upon  the  traffic 
engineer's  ability  to  control  the  traffic  within  the 
Lincoln  Tunnel,  some  of  the  theories  of  traffic  flow  will 
be  examined  in  the  next  section. 

1.4  General  Theories  of  Traffic  Flow 

Many  authors  have  studied  traffic  flow  in  detail  and 
only  a  brief  review  will  be  presented  here.   One  relationship 
basic  to  the  study  of  traffic  flow  is  response  at  time  (t+T) 
equals  sensitivity  multiplied  by  stimulus  at  time  (t).*   ' 
Response  is  usually  taken  as  acceleration  since  that  is  the 
variable  under  the  driver's  control.   The  stimulus  is  gener- 
ally represented  as  the  difference  in  the  velocities  of  two 


vehicles.   Thus,  £— «-  (X  (t+T))=  sensitivity  x  (4-X^(t)  minus 

.J.Z   n  at  n 

^— X  ^,(t))  where  X  represents  distance  traveled.   For 
dt  n+L 

mathematical  simplicity, sensitivity  was  originally  taken  as 
a  constant.   Later,  however,  more  sophisticated  models  of 


sensitivity  have  been  developed 


(12) 


A  fundamental  diagram  of  traffic  flow  has  been 
developed'*-'*)  and  is  shown  in  Fig.  2.   It  is  stated  that 


Density  (Vehicles/distance) 


>k 


FIGURE  2. 
Fundamental  Diagram  of  Traffic  Flow 

such  a  diagram  is  valid  only  for  a  particular  section  of 

roadway  at  a  particular  time.   Authors  have  placed  different 

interpretations  on  this  fundamental  diagram.   Martin^19' 

feels  that  it  is  valid  only  in  equilibrium  or  steady  state 

conditions  and  as  such  represents  the  desires  of  the  drivers. 

Duckstein  believes  that  Fig.  2  is  the  center  of  a  band  of 


6 
possible  flow-density  combinations.   Many  experimenters  feel 
that  traffic  density  (k)  is  the  best  indication  of  the 
actual  traffic  conditions^*  »21'.   Average  velocity  is  not  a 
good  indication  because  small  speed  changes  can  represent 
large  changes  in  flow. ^   '   Greenberg  has  compared  traffic 


flow  to  that  of  a  continuous  fluid. 


(13) 


Shock  waves  appear  to  be  a  major  cause  of  low  flow 
rates.   •  '   It  has  been  shown  that  a  small  increase  in  peak 
flow  rates  results  in  dramatic  reductions  in  traffic  conges- 
tton.<10> 

Authors  feel  that  density  can  be  used  to  predict 
the  formation  of  a  shock  wave^'  and  that  high  output  speeds 
indicate  the  absence  of  a  shock  wave.   Officials  of  the 
Tunnels  and  Bridges  Department  of  the  Port  of  New  York 
Authority  have  demonstrated  that  there  exists  a  critical 
output  velocity  and  that  traffic  densities  can  be  used  to 
foresee  sub-critical  velocities. dO)   Finally,  it  has  been 
shown  that  inherent  instability  exists  in  non -homogeneous 
traffic  flows.   '   The  velocity  at  which  traffic  actually 
travels  is  a  compromise  between  the  desired  velocities  of 
various  components  of  a  non-homogeneous  flow.   Thus  the 
components  which  are  required  to  travel  below  their  optimum 
velocity  can  react  to  a  sudden  decrease  in  the  speed  of 
leading  vehicles  in  such  a  way  as  to  create  a  shock  wave. 


7 
The  Port  of  New  York  Authority  has  been  and  is 
presently  attempting  to  induce  small,  frequent  gaps  in  the 
traffic  flow.   These  small  spaces  do  not  decay  as  rapidly 
as  larger,  less  frequently  spaced  gaps.   Only  a  small  degree 
of  success  has  been  realized  in  implementing  these  goals. 
The  south  tube  of  the  Lincoln  Tunnel  is  equipped  with  a 
closed  circuit  television  monitoring  system,  but  as  of  this 
writing,  the  tunnel  officials  only  react  to  but  do  not  pre- 
vent the  formation  of  shock  waves. 

The  traffic  sinks  in  New  York  City  and  the  Lincoln 
Tunnel  have  been  discussed  briefly.   This  is  necessary  for 
an  understanding  of  the  objective  of  the  thesis  which  is 
presented  in  the  next  section. 

1.5.   Objectives  of  the  Thesis 

Once  a  vehicle  departs  the  tunnel,  his  problems  have 
just  begun,   A  driver  will  find  himself  fighting  traffic 
all  the  way  to  his  particular  traffic  sink  and  then  fighting 
for  a  position  in  the  sink.   It  is  felt  that  by  requiring 
a  vehicle  to  wait  in  New  Jersey  before  entering  the  tunnel 
it  may  be  possible  to  allow  his  traffic  sink  to  depopulate 
sufficiently  so  that  the  total  time  required  for  him  to  be 
serviced  in  New  York,  is  actually  reduced.   The  basic  idea  is 
to  remove  this  driver  from  the  congested  area  long  enough 
for  some  vehicles  in  the  area  to  depart. 
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Since  it  is  almost  an  impossible  task  to  determine 
the  destinations  of  cars  and  small  vehicles,  it  was  decided 
to  apply  this  principle  to  large  trucks  and  tractor  trailers 
which  supply  and  remove  goods  from  the  C.B.D.  of  New  York 
City. 

The  general  model  is  presented  in  Chapter  2.   Two 
basic  problem  areas  arise:   (1)   sampling  the  sinks  to 
estimate  various  parameters  and  (2)  determining  if  the 
tunnel  supplies  a  large  enough  percentage  of  the  sink 
traffic  to  make  effective  control  of  the  sink  conditions 
possible.   These  problems  are  treated  separately  in  Chapters 
3  and  4. 


CHAPTER  2 
MODELING  THE  SYSTEM 

2.1.   General  Model 

As  a  vehicle  approaches  the  Lincoln  Tunnel  In  New 
Jersey,  It  enters  a  large  and  complex  traffic  system.   This 
system  consists  of  the  waiting  areas  in  New  Jersey,  the 
tunnel,  the  network  of  roads  in  New  York  leading  to  the 
destination,  and  the  traffic  sink  itself.    Upon  being  ser- 
viced and  departing  from  the  traffic  sink,  a  vehicle  is  con- 
sidered to  have  exited  the  system.   The  Lincoln  Tunnel  may 
be  thought  of  as  an  entrance  way  leading  to  24  separate  but 
not  wholly  independent  sub-systems  or  traffic  sinks.   As  a 
vehicle  leaves  its  sink   it  may  travel  into  another  sink 
and  thereby  become  part  of  the  uncontrollable  input  into 
another  sub-system 

Total  time  in  the  system  (T)  consists  of  six  compon- 

ents  (t.)  such  that  .E,t.  =  T  where: 
1  isL  i 

t  =  waiting  time  in  New  Jersey, 

t2  =  travel  time  in  the  tunnel, 

t~  ■  travel  time  from  tunnel  exit  to  desired  traffic 

sink, 
t^  =  waiting  time  in  sink, 
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tc  *  service  time  in  sink, 

t£  ■  time  required  to  exit  sink  after  being  serviced. 

The  time  required  for  a  vehicle  to  return  to  its 
original  destination  is  assumed  constant  and  is  not  con- 
sidered in  the  model.   All  times  are  measured  in  minutes. 
This  model  is  valid  for  any  time  period  during  the  day,  how- 
ever, here  it  will  be  applied  to  the  morning  rush  hour 

traffic  between  the  times  t      and  t   J#   See  Fig.  &. 

start      end 

All  times  are  assumed  to  be  between  these  limits 

2.2.   Components  of  the  Model 

The  individual  components  of  T  will  now  be  investi- 
gated separately. 

2.2.1.  Waiting  Time  in  New  Jersey 

Waiting  time  in  New  Jersey  (t  ) is  the  independent 
variable.   T  will  be  minimized  by  controlling  t..   An  upper 
bound  of  25  minutes  has  been  chosen  for  this  variable. 

2.2.2.  Travel  Time  in  the  Tunnel 

It  was  felt  originally  that  in  order  to  predict  t 
satisfactorily  it  would  be  necessary  for  the  model  to  predict 
the  status  of  each  vehicle  in  the  tunnel  at  every  second. 
This  involved  estimating  distance  traveled,  instantaneous 
velocity  and  acceleration  and  headway  for  every  vehicle.   It 
soon  became  obvious  that  this  approach  would  be  of  little 
practical  value  because  of  the  large  amount  of  computer 
memory  space  and  computing  time  required  to  predict  t2» 


LL 


To  avoid  these  problems,  it  was  decided  to  utilize 
the  flow-density  data  available  for  the  tunnel.  '     A  quad- 
ratic curve  was  fitted  to  this  data  by  the  method  of  least 
squares,  and  the  solution  of  the  normal  equations  was: 

9  a  a  ♦  bk  ♦  ck^  where 

9    ■   flow  measured   in  vehicles/minute, 

k  =   concentration  measured   in  vehicles/mile, 

a   =  0.44,    b  =   0.61,    and   c  =-0.0043. 
The  predicting   curve   is   shown   in   Fig.    3. 
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Flow  (e)  vs  Concentration  (k) 
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12 
It  Is  also  possible  to  observe  flow  as  a  function  of 
time  of  day.   This  data  was  obtained  and  the  flow  between 
the  times  of  t  t  t  and  t  ^  was  approximated  by  a  quadratic 
function.   This  curve  is  shown  in  Fig.  4,   The  quadratic 


e 

r*N 

a> 

/i 

3 

20- 

C 

•f4 

6 

1S 

CO 

<D 

V 

10-. 

«•* 

O   0) 

5- 

IM  > 

l"4  V^ 

t 1 r— i — i 1 1 1 1 1 1 1 1 

2    4 ,     6a  8    10  12    14   16   IS   20   22   24 
^tartVSnd  *ime(Hours)   <>t  Day 


t* 


FIGURE  4. 
Flow  (e)  vs  Time  of  Day  (t) 


equation  is  0  =  d  t2  +  pt  +  q  with  d  =  -9.2,  p  =  127.8  and 
q  «  -241. 

If  it  is  assumed  that  no  shock  waves  will  appear  in 
the  traffic  flow,  then  9  remains  fairly  constant  for  any 
given  time  of  transit  through  the  tunnel.   For  any  time  of 

the  day, 

2 

ck  +  bk  +  a  •  e  s  0, 


k  ■ 


-b+vH^. 


4c(a-0) 


2c 
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k  =  *>6^.61^4(>.  0043)  (.44-9) 
2(-.0043) 

i'or  t        t  -  t   ..  6  will  always  be  greater  than 
start  -   -  end»  ' 

.44  but  the  expression  under  the  radical  will  always  be 

positive.   This  can  be  seen  by  substituting  9     for  0. 

J  "  max 

Since  the  model  predicts  traffic  concentration  for  rush 
hour  traffic,  it  seems  reasonable  to  assume  that  the  actual 
density  is  the  larger  of  the  two  densities  corresponding  to 
a  specific  value  of  G.   Thus, 


k  a  ,6lW.6^-4(-,Q043)(t44-9) 
2(.0043) 

The  average  velocity  through  the  tunnel  (v)  is  the  flow 
divided  by  the  density.  Thus,  velocity  in  miles/minute 
may  be  written  as; 

v  =  e/k,     

v  =*  29c/(-b^b2-4c(a-9)).  (1) 

As  expected,  Eq.  1  simply  says  that  on  the  heavy  concentra- 
tion half  of  Fig.  4,  velocity  increases  then  decreases  with 
e. 

It  was  shown  that  9  ■  dt*  +  pt  +  q  for  tgt   t  <  t  < 
t   .    This  expression  is  substituted  into  Eq,  1  to  obtain 
velocity  as  a  function  of  time: 

v  ■  2c(dt2+pt+q)/(-bWb2-4c(a-dtz-pt-q)). 
Since  this  velocity  is  considered  to  be  the  average  velocity 
through  the  tunnel,  the  expected  time  required  for  a  passage 
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through  the  tunnel  is  found  by  dividing  the  tunnel  length 
(L  miles)  by  this  average  velocity.   Thus,  after  inserting 
the  values  of  the  parameters, 


t,  =  L( .  6l4^76p-4U,0043)  ( .  44+ 2. 36t2-33t+97)  )  (2) 

(-2.36tZ+33t-97)(2)(.0043) 

Since  t,  the  time  at  which  a  vehicle  enters  the  tunnel,  is 
the  time  of  arrival  at  the  New  Jersey  entrance  plus  the 
waiting  time  in  New  Jersey,  it  is  seen  that  t«  depends  on 
the  control  variable  t..   As  t  varies  between  tgt   t  and 
t   ,,  e  increases  and  then  decreases  (Fig.  4).   This  results 
in  the  unlikely  conclusions  that  t?  reaches  its  minimum 
value  when  t  «  t*t  or  when  6  is  at  its  maximum  value.  This 
is  shown  graphically  in  Fig.  5  in  which  t«  is  plotted  as 
a  function  of  t  as  t  varies  between  t      and  tQn(^, 

This  erroneous  result  was  cuased  by  the  two  follow- 
ing unrealistic  assumptions: 

1.  The  actual  concentration  is  always  the  larger 
of  the  two  concentrations  corresponding  to  a  specific  9, 

2.  No  shock  waves  will  be  induced  in  the  flow. 
Although  it  is  true  that  during  peak  traffic  the  flow  is  on 
the  heavy  concentration  side  of  Fig.  2,  it  is  not  reasonable 
to  assume  that  the  flow  started  on  that  side.   As  the  morning 
rush  hour  approaches,  flow  increases  along  the  left  hand  side 
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FIGURE  5. 
Time  in  Tunnel  (t2)  vs  Time  of  Day  (t) 

of  Fig.  3.   At  some  point,  the  flow  shifts  over  to  the  right 
hand  side  of  the  curve;  this  can  happen  in  one  of  two 
ways; 

1.  Flow  can  reach  its  maximum  from  the  left  and 
then  decrease  down  the  right  side  of  Fig.  3, 

2.  As  flow  is  approaching  its  maximum  value,  a 
shock  wave  can  shift  the  flow  to  the  heavy  concentration 
side  of  the  graph.   Thus  the  true  expression  for  tj   as  a 
function  of  t  is: 
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t    *  ljJj±ijd^Cr.M^QQ&ni^h^ .  ?r ,- W . 3 1+42  ip  t         <t<9t 

(-9.2t2+127.8t-421)(2)(.0043)  start-       1 


or 


t     =  L(  •  6I^  *  6l/-4(-.  0043)  ( .  44+9 ,  2t* 
(-9.2t2+127.8t-42l)(2)(. 


-127,8t+421))      ;<t<t 


0043) 


■  end 


i 
whore  t  represents   the  time  at  which  maximum  flow  is 

achieved   if  no    shock  wave  occurs   or  the  time  of  a   shock  wave. 

The   solid   lane   in   rig.    6   shows  how  t2  varies   with  t 

if  no   shcok  vrave   is   induced;    the  dotted   line  traces   t2   if  a 

t 

shock  wave  is  encountered  at  t.   In  either  case,  it  can  be 

concluded  the  t2  is  a  non-decreasing  function  of  t. 

If  Fig.  4  were  used  as  a  guide,  one  would  decide 
that  the  flow  will  never  be  on  the  light  concentration  side 
of  the  flow-density  curve  until  evening  rush  hour.   This 
obviously  is  not  necessarily  true.   An  acceleration  shock 
wave  can  occur  as  9  is  decreasing  along  the  right  side  of 
Fig.  3.   If  this  happens,  t2  will  begin  to  decrease.   Fig. 
7  shows  this  penomena  graphically. 

2.2.3  Travel  Time  to  the  Traffic  Sink 
The  author  had  hoped  to  study  t~  in  detail,  but 
because  of  time  limitations,  this  was  impossible.   Therefore 
t^  will  be  taken  as  a  constant  for  each  traffic  sink.  Refer- 
ences (2)  and  (23)  among  others,  discuss  some  of  the  prob- 
lems that  would  be  encountered  in  predicting  t^. 
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2. 2. A  Time  in  the  Traffic  Sink 

It  is  almost  impossible  to  know  exactly  what  will 
happen  to  a  vehicle  once  it  enters  a  traffic  sink.   Some 
vehicles  will  find  a  legitimate  servicing  station  immediate- 
ly, some  will  create  their  own  servicing  stations  by  double 
parking  or  by  parking  in  a  restricted  zone,  and  some 
vehicles  will  be  forced  to  wait  until  a  servicing  station 
becomes  available.   It  is  also  impossible  to  determine 
exactly  what  will  happen  to  a  vehicle  once  it  locates  and 
enters  a  servicing  station.   Some  vehicles  will  find  men 
available  to  unload  or  load  cargo  while  other  vehicles  will 
be  required  to  wait  for  servicers.  As  the  congestion  mounts, 
servicers  will  begin  to  interfere  with  each  other  and  will 
increase  service  time.   It  becomes  clear  that  the  standard 
queueinc  models  are  not  adequate  to  describe  the  situation. 

It  is  hypothesized  that  the  total  time  spent  in  a 
traffic  sink  (t.  +  tc  +  t,)  increases  with  the  number  of 
vehicles  in  the  sink  at  the  time  of  arrival  (N).   It  is 
further  hypothesized  that  this  increase  is  not  linear  but 
that  the  rate  of  change  of  the  total  time  in  a  sink  with 
respect  to  the  number  in  the  sink  also  increases  with  N. 

Letting  w  =  t,  +  tr  +  t6,  the  following  model  will 
be  hypothesized:   E(w)  =*  ae^1  ,  where  a  and  3  are  parameters 
which  must  be  estimated  for  each  traffic  sink.   The  possi- 
bility of  using  other  models  for  w  is  discussed  in  Chapter  3. 
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for  the  moment,  it  is  assumed  that  t~  =  t^  =  0  so 
that  vehicLes  enter  their  sinks  directly  and  instantaneously 
from  the  New  Jersey  waiting  areas.   Total  time  in  the  sink 
will  now  be  studied  as  a  function  of  t, . 

If  it  is  assumed  that  all  entries  into  a  sink  occur 
instantaneously  at  the  beginning  of  each  minute,  then  a 
consequence  of  the  model  is  that  all  the  vehicles  that  enter 
a  particular  traffic  sink  during  the  same  minute  have  the 
same  expected  tine  in  the  sink  ^(w).   Trucks  which  enter  a 
traffic  sink  during  the  same  minute  are  said  to  comprise  a 
block  of  vehicles.   At  time  t,  a  given  sink  may  be  represent- 
ed as  in  Fig.  8.   Here  N.  represents  the  number  of  vehicles 
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FIGURE  8 
Conditions  within  a  Traffic 
Sink. 


expected  to  depart  at  time  t..   The  expected  number  of 

entries  per  minute  from  all  sources  other  than  the  tunnel 

is  denoted  by  A. .    It  is  assumed  that  \    remains  constant 

for  tgtart< t< tend  and  that  the  first  in  first  out  rule  is 

valid  for  the  traffic  sinks.   Letting  N  ■  £  N<  the 

i=a  x 

following  relationships  can  be  devloped: 


if   t     =0,    E(v)    =  ae 
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B(N+X+1) 


tL  -   1,    E(w)    =  ae^+2X+1> 


,      /    \  0(N+(t,+l)X+l),        „     ^  *. 

and   in  general   l(w)   ■  ae  1  for  t^  <   t     -   t. 

The  sum  t,    +  E(w)    =  t,+ae'  ^      *    i      '        '    reaches    its  minimum 

value  when    t^   =   0. 

For   t     -   t  <  t.    <  t.    •   t     The  general   expressions   for 
a  "~*      *  d 

expected   time   in   the  sink   is   E(w)    =  ae^N"Na+^lH^+1\      The 
sum  of  t,    and  w  is  minimized  when   t,    »  ta  •   t.      Similarly 
for  t,    •   t  <  t|   <   t.   •  t   the   total   time   in   the  system   is 
t^e^N-^^b-'-NH-Ct^UX+l)      and   ts  minlini2ed  when   tL  - 


*1  '   C' 


It  is  also  seen  that  the  value  of  t,  which  minimizes 


T  is  such  that  a  vehicle  will  arrive  at  the  traffic  sink  at 
the  expected  time  of  departure  of  a  block  of  vehicles  and 
never  between  the  expected  times  of  departues  of  two  blocks. 
Unfortunately  to  ^  0  and  t~  j*  0  and  it  is  not  poss- 
ible to  control  t.  so  that  a  vehicle  will  enter  a  traffic  sink 
just  as  a  block  of  vehicles  is  departing.   Considering  that 
vehicles  must  go  through  the  tunnel  before  entering  their 
desired  sinks,  it  is  clear  that  the  present  expression   for 
w  is  not  adequate.  When  a  vehicle  arrives  at  the  tunnel,  the 
control  personnel  must  predict  the  conditions  at  his  desired 
traffic  sink  for  the  coming  25  minutes.   Then  a  decision  can 
be  made  as  to  how  long  to  delay  the  vehicle  in  New  Jersey. 
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This  prediction  must  consider  the  vehicLes  bound  for  this 

specific  sink  which  have  already  entered  the  tunnel  but  have 

not  yet  arrived  at  the  sink. 

Let  it  be  the  number  of  these  vehicles  which  have 
o 

passed  the  control  personnel  but  have  not  yet  arrived  at 

the  sink.   Then  the  expression  for  L'(w)  must  include  a  term 

f(M  ,ttt  )  which  predicts  the  arrival  of  these  vehicles, 
o   '  1 

Let  g(N,t,t.)  predict  the  dissipation  of  vehicles  from  the 
sink.   The  expression  for  E(w)  is  now 

£(w)  *  aexpO(No-g(N,t,t1)+f(Mo,t,tL)+(t1+l)X-H)). 

2.3.      Minimizing;   Total   Time^in   the  System 

Combining  all  parts  of   the  model  yields:      T(t,t^)=t^+ 


L(.61-rv/.6l^-4(.. 0043)  (. 44+9.  2(t+tL)2-127.8(t+T   )+421)) 
(.9.2(t+t   )    +127.8(t+t1). 421(2) (.0043) 

+   constant   +• 

0(N-s(Nft,t1)-ff(^o,t,t1)+(t1+l)X+l)  (3) 


where, 


016 


t   =   time  vehicle  arrives   at   the  tunnel 

t 


r  m\*l>   t.t«rtf£<e«i>   <  « 
(-1,   I  <  (t+tl)   <  tend. 

Chapter  3  discusses  ways  of  estimating  a  and  3. 

Assuming  for  now  that  this  has  been  done,  there  are  at  least 

two  ways  of  finding  the  value  of  t.  which  minimizes  T. 
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1.  3y  assuming  the  form  of  f(M  t,t  )and  gCN.t,^); 
it  is  then  possible  to  solve  Eq.  3  numerically  for  the  opti- 
mum value  of  t  .   An  example  of  this  is  done  in  Chapter  4. 

2.  By  computer  simulation;  the  case  for  this  method 
is  strong.^  »20fl°)   xhe  tunnel- sink  system  is  extremely 
complex  and  it  is  impossible  to  describe  mathematically 
exactly  what  wilL  happen  to  a  vehicle  which  enters  the 
system.   With  simulation  it  is  not  necessary  to  guess  f  and 
g,  nor  is  it  necessary  to  assume  first  in,  first  out.   As  a 
vehicle  enters  the  tunnel,  its  expected  time  in  the  tunnel 
is  calculated.   When  this  predicted  time  is  reached,  the 
vehicle  is  simply  transferred  out  of  the  tunnel  into  the  net- 
work of  roads  leading  to  the  sink  and  its  time  of  arrival  at 
the  sink  is  predicted.   It  has  been  assumed  that  a  decision 
is  made  for  only  one  vehicle  at  a  time.   With  simulation  it 
is  possible  to  determine  the  effect  of  letting  any  number  of 
vehicles  through  to  the  same  sink.   Time  of  arrival  at  the 
tunnel  (t)  is  really  a  random  variable.   With  simulation, 
the  modeler  can  study  the  effect  of  variations  in  t  as  well 
as  variations  in  the  initial  conditions  in  the  tunnel  and  in 
the  sinks.   Such  a  program  has  been  written.   The  program 
listing  is  in  Appendix  A  and  the  logic  flow  chart  is  in 
Appendix  B. 
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Before  the  experimenter  can  apply  the  theories  of 
this  chapter  to  a  particular  system,  he  must  be  satisfied 
that  he  can  predict  the  expected  waiting  time  in  the  traffic 
sinks.   This  problem  is  discussed  in  the  next  chapter. 


CHAPTER  3 
DRAWING   SAMPLES   FROM  THE  TRAFFIC  SINKS 

3.1.      ..^ppnential  i-iodel 

As  was  discussed  In  Chapter  2,  it  is  assumed  that 
the  sum  of  waiting  time  in  a  traffic  sink,  servicing  time 
and  exit  time  (w)  is  equal  to  ae  .  For  any  specific 
traffic  sink,  the  experimenter  must  be  able  to  estimate  a 
and  g  and  obtain  some  idea  of  how  well  his  model  represents 
the  conditions  of  the  physical  sink.  The  purpose  of  this 
chapter,  then,  is  to  present  a  sampling  technique  by  which 
these  aims  can  be  achieved. 

To  be  specific,  the  model  for  total  time  in  the  sink 
is: 

3N    r  9s 

w.  =  ae   +  £ .  where  t  ^N(0,6  )  ,  i  =  1,  2,  •  •  •  m, 
and  m  «  the  number  of  observations. 
It  is  assumed  that: 

a.  deviations  are  normally  distributed, 

b.  deviations  have  a  common  variance  6%  and 

c.  deviations  are  independent. 

The  actual  graph  of  w  vs  N  depends,  of  course,  upon 
the  specific  values  of  a  and  0,  but  In  general, it  can  be 
represented  as  in  Fig.  9. 
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Time  in  Traffic  Sink  (w)  vs  Number  of  Vehicles 
in  Traffic  Sink  (N) 


Since  tf(0,62)  a  .1 —  exp(-V(w-ae0N)2»6"2) ,  the  like- 

iihood  function  is  L»(JL)m/2.-*cn*exp(-   2.  (w  »ae   i)  .6   ). 

2rr  6  i=l      i 

Since  N.    and  Wj    are  observable  values,    there  are   three  para- 

2 
meters  which  must  be  estimated:    a   and  £,    and  6    .      L   is  maxi- 
mized by   selecting  a  and  0    such  that: 


m 


i=l      i  i=l 


(4) 


and 


i»l     ii  1*1      i 


(5) 
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Solving  Eqs.  4  and  5  simultaneously  yields  the 
maximum  likelihood  estimates  of  a  and  3.   It  is  interesting 
to  note  that  these  are  identical  to  the  least-squares 
estimates  of  a  and  p.   The  sum  of  the  squared  deviations  of 

the  observed  values  of  w^  from  the  predicted  values  is 

m        RN   2 
S(a,g)  ■  Z      (w^ae0  i)  .    Upon  taking  the  partial  derivi- 

tivies  of  S  with  respect  of  a,  and  0  and  equating  them  to 

zero  to  minimize  S,  Eqs.  4  and  5  again  result. 

To  obtain  the  maximum  likelihood  estimate  of  62  one 

9 

simply  maximizes  the  partial  of  L  with  respect  of  6  .   The 
result  is  6   ■  £   (w^-aep  1)  /m  where  m  is  the  number  of 
(w^,N^)  two-tuples  observed  and  •  over  a  parameter  represents 
an  estimate  of  that  parameter.   It  is  not  possible  to  deter- 
mine whether  or  not  o  is  biased;  this  requires  a  knowledge 
of  the  bias  in  a  and  p.   However,  maximum  likelihood 
estimators  are,  in  general,  not  unbiased. 

One  advantage  of  maximum  likelihood  estimators  is 
that,  under  fairly  general  conditions,  it  can  be  shown  that 
they  are  efficient  estimators.^6'   From  this  it  follows 
that  maximum  likelikhood  estimators  are  also  consistent. 

3. 2.   Determination  of  Sample  Size 

At  this  point  it  would  b e  desirable  to  have  a  back- 
log of  samples  from  which  it  would  be  possible  to  obtain  an 
advance  estimate  of  experimental  error  (6  ).   From  this  the 
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experimenter  could  determine  the  sample  size  required  to 
satisfy  his  precision  requirements.   If  such  a  backlog  of 
samples  is  not  available;  it  will  be  necessary  to  obtain  an 
advance  estimate  of  6^  in  another  way. 

The  proposed  methods  of  drawing  samples  and  estimat- 
ing parameters  will  be  demonstrated  by  means  of  an  example. 
The  traffic  sink  selected  was  a  Circle  K  food  market  in 
Tucson,  Arizona.   It  is  not  hypothesized  that  the  Circle  K 
represents  a  New  York  traffic  sink;  however,  the  food  mar- 
ket does  possess  some  quantities  similar  to  those  of  a 
traffic  sink.  First,  most  people  enter  Circle  K  with 
similar  missions  --  to  buy  one  or  two  items.   In  ;Iew  York, 
trucks  entering  traffic  sinks  have  similar  missions  also. 
Second,  there  is  reason  to  believe  that  the  congestion 
created  by  the  presence  of  many  cars  in  a  Circle  K  park- 
ing lot  will  cause  the  rate  of  change  of  waiting  time  to 
increase  with  N. 

For  every  entering  vehicle,  it  is  necessary  to 
determine  w.  and  N. .   A  convenience  sample  of  size  m«8  was 
drawn  and  the  results  are  shown  in  Table  1. 

A  modified  steepest  descent  method  using  an  itera- 
tive search  was  used  to  solve  the  nonlinear  equations  1  and 
2.   The  following  results  were  obtained: 

a  »  0.726,  p  =  0.64  and  52  ■  0.190. 


TABLE   1.      Results   of  Pilot   Sample 
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i 

r     ~   Ni 

Wj    (minutes) 

1 
2 

C 

1 

iSi  Ni=8 

1 

2 

3 

2 

3,    2,    2,    3 

4 

3 

5,    5 

iiow,   using   this    estimate  of   experimental   error,    it 
is  possible  to  determine  the   sample  size  necessary  to   satis- 
fy  the  precision   requirements.      These  requirements  may  be 
worded   in   several  ways.      A   few  examples  will   be  given. 

3.2.1.      Requirement    I: 

Determine  the   sample   size   (m) ,    such  that   S(w)    can  be 
estimated  within  \  minute  99.5%  of  the  time. 

Since   (T  QQc;  ,)   £  =  d  =  %  minute, 


m 


.995,  m-l/y= 
<t.995.m.l)2«2>   ,„ 


. 995, m-1 


)Z(0.76). 


It  is  necessary  to  iterate  on  m  until  the  equality 
is  satisfied.   It  was  found  that  with  m=8,  (t  gg5   7)2(,076) 
equals  9.75.   with  m=9,  (t  995  g)(0.76)  ■  8.51.   Therefore, 
for  these  precision  requirements  a  sample  size  of  9  is 
necessary. 
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3.2.2  Requ ir enen  t  II: 

Determine  the  sample  size  such  that  C  a  99%  confi- 
dence limits  on  2(w)  will  have,  with  probability  p=.99,  a 
length  less  than  or  equal  to  one  minute. 

Here  we  use  the  relationship  developed  by  Harris, 
Harvitz,  and  Mood:^*5' 
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/.  1 2 


where 


v+1  ■  pilot  sample  size  =  7,  and 
d  »  1  minute. 

Hence,  1.32  „  ^.P^Lm-I^Ol."-!.^  . 

m 

Again   it   is  necessary  to    iterate  an  m  until   the  equality   is 
satisfied.       In   this   example,    a   sample  size  of   35  will   fulfill 
the  precision   requirements. 

3.2.3.      Requirement    III: 

Determine  the   sample  size  required   for  the  probabil- 
ity of   a   90%  confidence   interval  on   £(w)   having   a   length 
d<6/5  to  be  equal   to   95%. 

Letting  w.    be  the  estimate  of   ^(w^)    it   follows   that: 

P<  ?!•£<*].))    <   8ttg5,ni-l  =  Q#90  and 

»/m 


P(   w^cCwj^))    <  6/10)    =  0.95. 


Hence, 


st    Q? 

■^  <  6/10. 
m 


30 
To  find  the  minimum  value  of  m  that  satisfies  the  given 

requirements,  set    .95tm-l  =  6/10. 

Vm 
Squaring  both  sides  and  simp lying,  Eq.  6  is  developed: 

ns^   m- 

~^   l00cZ0.95,m-l-  (6) 

It  is  now  observed  that  SSJ.^  )C*" 

6  2     m-1 

For  values  of  v  greater  than  30,  it  is  possible  to 
use  the  fact  that  C/2%2  -  V2V-1)2.  is  distributed 
approximately  as  N(0,1)  to  obtein  the  following  approximation, 

Now,  incorporating  this  result  into  ;:q.  6  the  result  is, 

0C2. 95  v  "  %(Z.95  +V*?rT)2=  mi_^ 

100t  .95,-1   • 

A.95,n,.l      -95  lOOt^js^!   , 

hence 

m  2=  50t2.95,n>.l  <2.95  +  <&*>*' 
For  values  of  m  greater  than  30,  the  approximation 

t.95fm-lrvZ.95  is  VaLid' 

Therefore,  m  ■  50Z"  qs(S  Ae,  -K/2m«3)^  and  because  a  negative 

'  •     .  95 

sample  size  is  meaningless, 

V3o  z>95  (z>95  +  .ySnTS), 

«2   -    2(19. I5)m  +   (19. 15)2   =   (ll.62)2(2m-3) , 


m 


m 


ra2   -    310. 3m  +  775  +  0. 
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Therefore: 


m 


m   +310.3W310.32.4(1)(775) 


31 


31C.3  ±   303 


m  =  305  or  4.   Upon  substituting  these  two  sample  sizes 
into  the  precision  requirements  it  is  found  that  m  =  4  is 
too  small  and  m  ■  305  is  the  only  acceptable  solution. 

It  is  readily  seen  that  Requirement  III  is  by  far 

the  most  restrictive;  it  does  not  require  an  advance  esti- 

2 
mate  of  5  ,   It  is  felt,  however,  that  few  people  would  word 

their  precision  requirements  in  this  manner.  Therefore,  for 

this  example,  Requirement  II  was  selected  and  a  systematic 

sample  with  a  random  start  of  size  m  =  35  was  drawn.   Ths 

results  are  summarized  in  Table  2. 


TABLi;  2.   Results  of  Sample 


I 
1 

2 
3 
4 

5 

°i 

"i 

1.5", 

Wij 

Hi 

,sr  w. 
j=i 

10.0 

8 
8 

0 

1.0, 

0.5, 

1.0,    2.5,    1.0, 

1.0, 

1.5 

1 

2.0, 

2.0, 

1.5. 

3.0,    1.5,    2.0, 

i.o, 

2.0 

15.0 
21.0 

9 
7 

2 

2.0, 

2.5, 

2.5, 

2.0,    2.0,    3.0, 

2.0, 

2.0, 

3.0 

3 

3.5, 

3.5, 

4.°. 

4.5,    4.5,    3.5, 

3.0 

26.5 

3 

4 

6.0, 

5.0, 

6.0 

17.0 
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3.3.      Es t ima t ion^o  f m P  a ramet er s 

It    is   convenitent  to  use  a  and  3    from  the  pilot 
sample  as    initial  guesses   in   the   iterative  technique  for 
solving   £qs.    4  and  5  for  the  final   estimates  of  a  and  3. 

The  results   obtained  were: 

1  •  »2  2 

a  »  0.591  minutes,   0    ■  0.590,    6      =  0.599  minutes    , 

i 

6    =  0.774  minutes, 

w.    -  0.591  minutes^         w^   =   1.066  minutes, 

w^  ■  1.922  minutes,    w.  =  3.468  minutes, 

w     »  6.256  minutes. 

The  experimenter   is   now  in  a   position   to   test   assumption  c 


in   section     3.1.    that   the  deviations  of  Wi    at   each  value  of 

ie  a1 
(1) 


o 
N^  have  a    common  variances  6    .      For  this   the  author  will  use 


Bartlett's   test  of  Homogeneity  of   Variance. 

For  this   test, 

k/(l-L)  «J  0A>1  where 

K  a   2.3026  ["(.£_    dfjlog   s2   -    (Zdf.    log   s   2) 
L    l"l        i  p  1  1    _ 

dfT  "  2df . 
3(k.l)      , 

k  ■   the  number  of  populations  of  w's, 

2     S3,   +  SS2  ♦  •••   +  SSi 

and  s   =  


Sdft 


In  this  example, 
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L  .  umgjm  .  0.08> 

s2   =  8.97/30  =  0.299, 
P 

log   s2  -  -0.524, 
P 

Zdf^   log    s^  =  -L6.15, 

K  ■  2.3026    (-15.72   +  16.15)    ■  0.99,   and 

X2   ■   .99/1.08   =  0.917. 

Since^Cr  =  9.48773  one  cannot   reject  HO:    62  =  62   ■ 

•  -'-'<"*  i.  ^ 

•••    =  6£  and   the  assumption   seems   to   be  valid. 

It   is   also  possible   to   calculate  95%  confidence 

limits   on  w. . 

t 
E(Wl)   =  w.    +  (t  95 >nl  >-*§*=  «i  ♦  df 

For  the  Circle  K   data, 

dL  -  0.52,  d2  =  0.52, 

d3  =0.48,  d     =  0.56,  and 

d5  -   1.31. 

Since  the  experimenter  is  using  the  same  estimate  of 
experimental  error  for  each  array,  the  length  of  the  confi- 
dence interval  depends  only  upon  n..   The  graph  of  w.  vs  I  . 
is  shown  in  Fig.  10  along  with  the  95%  confidence  at  each  w^. 
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w  =a+bN. 


x  indicates  observed  data 


*N 


12  3  4 

Number  of   Vehicles    in   Food  Market 


FIGURE   10. 
Time    in   Food  Market    (w)   vs  Number  of  Vehicles 
in   Food  Market    (w) 
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3.4.   iJva.luation  of  the  i-ipdel 

The  author  has  estimated  a,  3  and  6  for  the  traffic 
sink.   It  is  now  desired  to  examine  the  goodness  of  fit  of 

the  model,  to  the  physical  system.   experimental  error  (6)  is 

t 

always   an   indication   of  goodness  of   fit.      The  value  5    =C.774 

minutes  may  be  good,    it   is    impossible  to    say  how  good  how- 
ever   ,    since  there   is  nothing   to   compare   this  with.      As  was 
stated   earlier,    what   is   really  hypothesized  when   it   is 
assumed  that    i(w.)    =  cie^i   is   that  the  rate  of   change  of  w 
(as  well   as  w  itself)    increases  as  N   increases.      This   assump- 
tion  can  be  tested  by   comparing   the  model  with  a   linear 
regression  model. 

The   linear  model    is    E{wiJ    "d+gN^  +  D^+fj. 
where  w. .    is   the  jth  observation  from  the   i       population, 

D.    is   the  deviation  from  linearity  of  the  it3ri 

y-s  r,h 

population,   and  £ .    is   the  experimental   error  of  the  i 

population.      Again   it   is  assumed  that  the  experimental   error 

at   each  value  of  N.    is  normally  distributed  with  the  common 

variance  5^.      This  model  has   the  advantage  of   testing 

directly  the  hypothesis  that   the  deviations   from  linearity 

are  equal   to   zero.^      '      The  analysis  of  variance  for  this 

model,   using   the   same  data    is   given   in  Table  3. 

Since  E  (MS  departure  from  linearity)    =  62  -f  (k-2)"l+  L     N.D? 

and      E(MS   Residual)   =   52  we   can  test  HO:D.=D2='  *  ^Dj^O  with 
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TABLE  3.      Analysis   of   Variance   for  Linear  Model 


Source 

df 

SS 

MS 

Linear 
Regression 

1 

2lc                -  ? 
b  Z  m.(N.-N)    =52.26 
i«l    i      i 

52.26 

Departure   frOTt 
Linearity 

5-2   =   3 

k              _     m                -.9 
Z     n,(w  -w-b(N   -N)^=5.16 
i=l      i      i               i 

1.72 

Residual 

35-5=30 

k       m              —     9 
Z       Zl   (w.  .-v4)Z   =  ?.?7 
i=l    i=l        *-J      1 

0.299 

Total 

34 

k       ni    ,          «x2 

I        S1   (w.    -w)Z   =  66.39 

i=l   j  =  l        J-J 

55.54 

the  ratio 

F  *  i.S  Departure  from  linearity 
MS  Residual 

The  results  are  N  ■  1.69,  w  »  2.56,  6  =  0.92,  3  =  0.97, 
w.  =  1.25,  w~  *  1.87,  Wo  =  2.33,  w^  =  3.79  and  w<-  =  5.67. 

1.72 


F  = 


=  5.75,  and  F 


■  4.171.   Hence  HO  is 


0.299   »  —  *. 95,3, 30 

rejected,  and  the  author  concludes  that  the  linear  model  is 

not  adequate  for  predicting  E(w  ) .    It  has  not  been  shown 

that  ^(w. )  =>   ae°"i  is  the  best  possible  model  for  a  traffic 

sink  of  the  type  considered.   But  it  can  be  concluded  that 

it  is  better  than  the  linear  model.   The  linear  regression 

line  is  shown  in  Fig.  10. 
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If  the  exponential  model  had  not  given  a  better 
prediction  of  E(w.)  than  the  linear  model,  the  experimenter 
would  be  forced  to  change  the  model  of  the  traffic  sink. 

The  author  feels  the  polynomial  regression  would  do 
well  in  this  situation  since  the  experimenter  has  the 
ability  to  step  the  degree  of  polynomial  until  he  finds 
the  best  fit  of  the  polynomial  models. 

At  this  point,  it  can  be  concluded  that  the  expected 

waiting  time  of  vehicles  entering  the  traffic  sink  in 

6N 
question  can  be  represented  adequately  by  E(w^)  =  ae    . 


CHAPTER  4 

CONTROLLING  TRAFFIC  SINK  CONDITIONS 

4,1.   Introduction 

After  estimating  a  and  3  and  satisfying  himseLf  that 

he  can  predict  adequately  waiting  times  for  a  particular 

sink,  the  experimenter  still  must  answer  one  more  question 

before  the  philosophy  of  Chapter  2  can  be  applied  to  a 

tunnel- sink  system.   Let  y.   represent  the  expected  number 

of  arrivals/minute  at  sink  of  vehicles  which  have  traveled 

through  the  tunnel.   This  value  y.,  can  be  estimated  from 

the  origin  and  destination  data  available  for  the  tunnel  and 

must  be  the  average  over  a  certain  time  interval  within 

(t     .  ,  t   ,).   Again  letting  X,  equal  the  expected  number 
start   end  J 

of  arrivals  per  minute  from  all  sources  other  than  the  tunnel 
for  sink,  the  question  is  then,  does  the  tunnel  supply  a 
large  enough  proportion  y,/(Yj+^j),  of  the  expected  sink 
traffic  to  ensure  that  traffic  controls  at  the  tunnel  will 
be  reflected  in  the  sink? 

If  90%  of  the  traffic  in  the  sink  enters  via  the 
tunnel,  then  the  experimenter  knows  he  has  sufficient  con- 
trol over  the  sink  to  employ  the  method  of  Chapter  1.   iiow- 
ever,  if  only  1%  of  the  sink  traffic  comes  through  the 
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tunnel  the  experimenter  can  do  aLmost  anything  to  the  tunnel 
without  affecting  conditions  in  the  sink.   In  other  words, 
there  is  a  critical  value  k*   of  k   above  which  the  experi- 
menter does  not  have  sufficient  control  of  the  sink  condi- 
tions to  nerit  controlling  the  tunnel  traffic.   For  k>k* , 
it  is  expected  that  the  optimum  waiting  time  in  New  Jersey 
will  always  be  zero.   This  critical  value  of  k   must  be 
estimated  separately  for  every  sink  and  depends  upon  a,  3, 
the  expected  operating  conditions  of  the  sink,  and  the 
expected  rush  hour  conditions  in  the  tunnel. 

4.2,  Estimation  of  the  Critical  Value  of  k 

In  order  to  estimate  X*,  it  will  be  assumed  that 
f(M0,t,t.)t  the  function  which  predicts  the  arrival  of 
vehicles  which  have  left  the  tunnel  but  have  not  yet  arriv- 
ed at  the  sink,  and  g(ri,t,t^),  the  function  which  predicts 
the  dissipation  of  vehicles  initially  in  the  sink,  are 
known.   ^gain  let  M  be  the  number  of  vehicles  which  have 
entered  the  tunnel  but  have  not  yet  arrived  at  the  sink  and 

let  N  be  the  number  of  vehicles  presently  in  the  sink, 
o 

A  reasonable  first  approximation  for  f  and  g  is  to  assume 

that  they  decrease  linearly  from  their  initial  values.   If 

H   is  depleted  linearly  and  reaches  zero  in  U0  minutes  and 
o  ^ 

NQ  decreases  linearly  to  zero  in  U.  minutes  then: 
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~(^tt,tL)   -   (N0-lNOtl/Ul)kl         WherQ  Rl   =/ 


i,t1<u1 


o»Vui 


and 


f(M  tttl)  =  (M0-H0tl/U2)(T/H)k2 


where  "0  =  expiected  flow  for  the  next  U^  min.  , 

H  =  expected  no.  of  vehicles  in  the  tunnel  for  next 

U0  min. ,  and 


o,vu2  . 


Consider  the  second  equation.   The  term  (^0-;'0cj/U<?)  simply 

estimates  the  number  of  vehicles  bound  for  the  sink  which 

have  not  yet  arrived.   If  t,  ■  0,  this  term  is  M  and  if  t.=* 

U«  this  term  is  zero.   Since  t3  is  assumed  to  be  a  constant, 

vehicles  enter  the  traffic  sink  with  the  same  flow  with 

which  they  exited  the  tunnel.   If  all  *\  vehicles  were  ori  - 

o 

inally  in  the  tunnel,  then  (K  »M  t,/U„)/M  indicates  the  per- 

o  o  i  L 

centage  of  vehicles  in  the  tunnel  with  this  sink  as  their 
destination.   Since  u  vehicles/minute  exit  the  tunnel,  then 
(i:o-ii0t^/U2)(~/K)  vehicles  exit  the  tunnel  each  minute  for 
this  particular  traffic  sink,  and  t-  minutes  later  they 
enter  the  sink.   Then,  in  general, 

2  1  J 
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The  term  to  is  Introduced  to  account  for  the  time  delay 


between  exiting  the  tunnel  and  entering  the  traffic  sink. 
Notice  that  the  vehicles  which  enter  from  sources  other  than 
the  tunnel  are  never  removed  from  the  sink.   This  implies 
that,  besides  the  first- in,  first  out  assumption,  it  must  be 
assumed  that  these  vehicles  will  not  depart  from  the  sink 
before  t^  reaches  its  maximum  values.   Although  these 
assumptions  seem  reasonable,  they  are  not  strictly  necessary 
as  will  be  shown  later. 

4.3.   ^xample  of  estimation  Procedures 

To  demonstrate  the  method  of  determining  X*,  it  is 
hypothesized  that  a  scaled  dovn  version  of  the  Lincoln 
Tunnel  supplies  the  Circle  K  with  vehicles.   The  Circle  K  is 
used  as  a,  0  and  5   have  already  been  estimated  for  it.   It 
is  assumed  that  t3  =  1  minute  and  the  upper  limit  on  t^  is 
five  minutes.   Fig.  11  shows  the  plot  of  £(w)  vs  tj^  for 
various  values  of  X.   Time  in  the  tunnel  (t„)  is  simply  a 
scaled  down  version  of  Fig.  6.  without  shock  waves.   ..he 
values  of  t,  which  minimize  T  are  shown  in  Table  4,   For 
X   >  1.4,  the  optimal  waiting  time  (t.)  is  always  zero.   This 
is  because  for  X>1.4  the  experimenter  cannot  exert  suffi- 
cient control  over  the  traffic  sink  to  make  waiting  worth- 
while.  Uncontrollable  sources  overshadow  all  the 
experimenters  efforts  to  reduce  the  number  in  the  traffic 


Waiting  Time   in  New  Jersey 
(Minutes) 

figure  11 

Time  in  Traffic  Sink  (w)  vs  Waiting  Time  in 

New  Jersey  (t.) 


42 


X  =  1.5 


X=1.4 


sink.   Again  using  Fig.  6  as  a  guide,  a  similar  analysis 
was  made  assuming  that  a  shock  wave  occurred  in  the  tunnel. 
The  results  were  the  same  and  indicated  X*  -  1.4. 

In  order  to  circumvent  the  f irst-in-f irst-out  assump- 
tion  and  the  assumption  that  no  vehicle  which  enters  from  a 
source  other  than  the  tunnel  will  exit  the  sink  before 
t|  »  t,    ,  the  method  of  depleting  N  must  be  modified.   It 
can  be  assumed  that  a  constant  percentage  of  the  vehicles 
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TABLi^  4.   Linear  Depletion  of  f  and  g 


0 

.i 

o 

•  *— 

.3 

.4 

.5 

.6 

.7 

.8 

.9 

1.0 

1.1 

1.2 

1.3 

1.4 

1.5 


tj^  optimum 

T  minimum   (minutes) 

3 

3.9 

3 

9.4 

3 

9.9 

3 

10.6 

4 

11.2 

4 

11.9 

4 

12.9 

4 

14.1 

5 

15.4 

5 

17.1 

5 

19.4 

5 

22.5 

5 

26.6 

2 

30.3 

0 

31.0 

0 

31.0 

which  are  in  the  sinlc  at  the  beginning  of  one  minute  will 
have  edited  by  the  end  of  that  minute.   At  the  beginning  of 
each  minute,  vehicles  enter  the  traffic  sinlc  from  uncon- 
trollable sources  and  possibly  from  the  tunnel  and  combine 
with  vehicles  already  in  the  sink.   By  the  end  of  this 
minute  a  certain  percent  of  the  vehicles  have  departed. 
Besides  the  previously  mentioned  advantages,  this  approach 
does  not  stipulate  the  manner  in  which  vehicles  leave  dur- 
ing a  riven  minute.   This  method  does  require  an  updating  of 
N  at  the  end  of  each  minute  before  vehicles  enter  for  the 
next  minute,  but  still  falls  under  the  category  of  assuming 
the  form  of  f  and  g. 
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Fig  12  shows  w  as  a  function  of  t^  if  no  shod',  waves 
are  encountered  in  the  tunnel  and  if  the  constant  percent 


\*C.6 


Waiting  Time   in  New  Jersey 
(Minutes) 

FIGURE    12. 

Time  in  Traffic  Sink  (w)  vs  Waiting  Time  in 
New  Jersey  (t,) 

depletion  method  is  used  to  update  N.   ?or  k   >  0.6  it  is 
seen  (Table  5)  that  the  experimenter  does  not  have  suffi- 
cient control  over  the  sink  traffic  to  make  waiting  before 
entering  thd  tunnel  profitable.   Similar  results  were 
observed  with  a  shock  wave  in  the  tunnel.   T  minimum  was 
increased  but  X*  remained  0.6. 

The  experimenter  must  now  decide,  from  his  know- 
ledge of  the  traffic  sink,  which  method  to  employ  in 
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TA3 

LriCi       -* « 

Constant  Percentage  Depletion 

of   f  and  g 

X 

el 

optimum 

I  minimum   (minutes) 

0 

3 

8.3 

.1 

3 

8.9 

.2 

4 

9.3 

.3 

4 

10.6 

.4 

4 

11.8 

.5 

5 

13.7 

.6 

0 

15.9 

.7 

0 

15.9 

estimating  X*.   For  example,  the  author  feels  that  the  first 
in,  first  out  assumption  is  not  unreasonable  for  the  Circle 
K  sampled.   However,  he  does  not  believe  that  the  first 
method  can  be  realistically  applied.   While  making  a  deci- 
sion for  one  vehicle,  some  of  the  contontrollable  vehicles 


will  enter  and  exit  while  0  <  t,  <  t^ 


max. 


Although  this 


will  not  be  true  if  w  is  expected  to  be  larger  than  t.     , 
the  author  feels  that  for  the  Circle  K,  X*  =  0.6  is  a  more- 
reasonable  limit  than  X*  =  1.4. 

It  is  not  difficult  to  estimate  X.   While  sampling  to 
estimate  a  and  g ,  the  easiest  method  is  to  record  the 

initial  value  of  N  and  then  mark  the  time  of  each  arrival 

t 
and  departure.   For  this  X  can  be  determined.   Precision 

requirements  on  X  may  necessitate  drawing  a  larger  sample 
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than  is  indicated  in  Chapter  2,  however  no  theoretical 

problems  are  involved, 

i 

If  X  <  A.K,  the  experimenter  can  feel  that  he  has 

enough  control  over  the  sink  traffic  to  be  able  to  reduce 
T  by  controlling  t,. 


CHAPTER  5 
SUMMARY  AND  CONCLUSIONS 

5,1.   Summary 

Ihe  hypothesized  model  for  predicting  total  time  in 
the  system  is  T  =  t,  +  t2  +  constant  +  ae^N(tl).  As  TL 
increases,  N  and  hence  E(w)  decrease.   The  objective  is 
to  find  the  value  of  t.  for  which  the  reduction  in  T 
attributal  to  descreases  in  E(w)  is  offset  by  the  increases 
in  T  caused  by  the  growth  of  t  and  t2(t.). 

Before  attempting  to  use  this  model  on  an  actual 
tunnel-sink  system,  two  basic  questions  must  be  answered. 

1.  Can  the  total  time  in  the  traffic  sink  be 
predicted  adequately?   To  answer  this  question,  the  experi- 
menter must  sample  the  sinks  to  estimate  the  parameters  a 
and  3,  using,  for  example,  proceedures  outlined  in  Chapter  3. 

2.  Does  the  tunnel,  over  which  the  experimenter  has 
control,  supply  a  large  enough  percentage  of  the  input  to 
the  sink  so  that  the  reduction  of  M  will  overcome  the 
increases  in  T  caused  by  t,  and  t2?   This  problem  is  dis- 
cussed in  Chapter  4  and  involves  estimating  A.  and  X*. 

For  the  tunnel,  it  was  found  that  t~  increases  with 
t.  as  in  Fig.  6  of  Chapter  2.   For  the  Circle  K  it  was  found 
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that  a  ■  0.59,  0  =  0.591,  \   =  0.47,  and  X*  =  0.6.   Thus,  if 
a  scaled  down  version  of  the  tunnel  supplied  vehicles  to 
the  Circle  Kt  the  optimum  t,  would  be  4  minutes  and  would 
result  in  an  expected  total  time  in  the  system  of  13.3 
minutes.   These  results  are  summarized  in  Fig.  13. 


tl*ft2+t3+w 


lmax 


Waiting  Time   in  Hew  Jersey   (Minutes) 


FIGURE    13. 

Total  Time  in  System  (T)  vs  Waiting  Time  in 
New  Jersey  (t,) 

Since  it  is  reasonable  to  believe  that  t~  will  actually 
decrease  as  t   increases  and  thus  help  decrease  T,  the 
assumption  t«  a  constant  is  a  conservative  one. 
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5.2.  Application  of  the  i'odel  to  the  Lincoln  runnel  System 

There  are  several  stumbling  blocks  which  must  be 
overcome  before  the  idea  can  be  applied  successfully  to  the 
Lincoln  Tunnel  -  New  York  City  system. 

1.  There  are  no  vehicle  waiting  areas  available  in 
New  Jersey. 

2.  As  of  this  writing,  the  authorities  had  not  yet 
set  up  a  traffic  control  in  the  tunnel  so  as  to  avoid  shock 
waves . 

3.  Because  of  the  sizes  of  the  traffic  sinks,  the 
number  of  vehicles  involved,  and  the  large  number  of  sink 
entrances  and  exits,  the  sampling  problems  in  New  York  are 
expected  to  be  very  difficult,  much  more  so  than  for  the 
example  developed  in  Chapter  3. 

Due  to  the  myriad  of  entrances  into  any  New  York 
sink  the  question  of  controllability  posed  in  Chapter  4  is 
very  real.   It  is  unlikely  that  the  Lincoln  Tunnel  supplies 
a  large  enough  percentage  of  the  traffic  into  the  Central 
Business  district  of  New  York  City  to  make  this  method  of 
control  applicable  to  all  24  traffic  sinks, 

5.3.  Extension  to  U.  S.  Naval  Ships 

The  theory  and  methods  discussed  in  Chapter  2  apply 
to  any  transportation  system  in  which  controlling  the 
vehicles  as  they  enter  will  affect  the  conditions  within 
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the  system  that  these  vehicles  will  encounter.   To  the 
author,  one  of  the  most  obvious  extensions  of  these  control 
procedures  is  to  fiaval  ships  entering  port  through  the  con- 
straints of  a  channel.   Many  of  the  problems  encountered  in 
New  York  are  alleviated  in  this  channel-harbor  system  and 
some  new  difficulties  arise,  but  the  basic  control  theories 
remain  the  same. 

The  traffic  sinks  within  the  harbor  are  the  deep 
draft  piers,  the  shallow  draft  piers,  fueling  piers, 
anchorages,  possibly  ammunition  loading  areas  and  possibly 
several  points  within  a  shipyard.   There  are  fewer  sinks 
than  in  New  York  and  they  are  more  distinct  and  easier  to 
define.   The  sampling  procedures  are  simpler  than  those  in 
New  York  and  probably  would  require  little  more  than  sampl- 
ing the  watch  officers'  logs  of  selected  ships. 

There  would  be  no  requirement  to  construct  waiting 
areas  at  the  channel  entrance,  and  the  problem  of  control- 
lability is  completely  removed.   The  Senior  Officer  Present 
Afloat  (SOPA)  has  absolute  control  over  all  vehicles  which 
would  enter  the  traffic  sinkB.   Thus  X  would  be  zero  and  X* 
non-existent.   3ecause  of  the  speeds  and  spacings  involved 
it  is  much  easier  to  predict  conditions  in  a  channel  than 
to  predict  conditions  in  a  tunnel. 

The  interactions  between  ships  attempting  to  be 
serviced  (to  moor)  in  the  same  sink  are  just  as  real  and  as 
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bothersome  as  the  delays  caused  by  trucks  trying  to  off  load 
cargo  in  the  same  area.   It  is  believed  that  because  of  the 
simplifications  inherent  in  a  channel-harbor  system,  the 
control  procedures  outlined  in  the  preceding  chapters  can 
be  applied  successfully  to  these  systems  at  this  time. 

5.4.   Conclusions 

The  model  of  Chapter  2  is  by  no  means  complete.   For 
example,  no  mention  has  been  made  of  the  Holland  Tunnel, 
another  facility  under  the  control  of  the  Port  of  New  York 
Authority. 

The  initial  study  of  a  large  problem,  must  by 
necessity,  involve  many  simplifying  assumptions.   These 
enable  the  experimenter  to  isolate  major  problem  areas  as 
was  done  in  Chapters  3  and  4. 

In  conclusion,  it  is  believed  that  the  idea  of 
reducing  traffic  delays  in  New  York  City  by  controlling 
vehicles  before  they  enter  the  tunnel  is  possible  and 
worthwhile.   Authority  officials  feel  that  it  will  be  years 
before  they  actually  begin  experimenting  with  the  ideas 
developed  in  this  thesi6.   But  the  problem  of  traffic  con- 
gestion ,  outlined  in  Chapters  1  and  2,  must  someday  be 
faced.   It  is  hoped  that  this  study  constitutes  a  first 
step. 


APPiH?DIX  A 
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00  720  J  =  1,5 

TIS7  =  TIS7  +  ADO 

WT  s  AA  (  I  A  )  *E*  PEF ( P ( I  A ) *T I  57 ) 

3M{  II  ,J)  =  BM<  I  I  •  J )  +  WT 

ADD  '-     ADD  +  1.0 

5=5+1.0 

TIS7  -    TIS7  -  ADD 

DO  050  II  =  1 ♦ 2 5 

DO  0  50  J  =  ] »5 

IBMt I  I »J)  -  BM{  II tJ 

DO  9960  I  =  1 »25 

FORMAT ( I  X »5I6 ) 

Pr;  INT  9961  *     i   IBM!  I  ♦ 


.5 


) 


J  =  1 


IBM (  I  I  ,1  )  ] 

1 1 1 » J ) 


DO  o52  II  =  1 »?' 
I  COL  (II)  --  1 

DO   052  J  =  2.5 

if  f i°"( r i « J)  - 

I  B«  ( 1 1  » 1  )  =  I P  M 

I  COL  (II)  =  J 

CONT INUE 

NCOL  •--  I  COL  (  1  ) 

NPCVt    =  1 

DO  954  T  T  =  2.25 

IF  (IBM(II.l)  -  IBMtltDJ  955.955.954 

IBM* 1 .1  J  =  IBM  (  I  I  ♦  1 ) 

NCOL  =  I  COL  (II) 

NPO'\i  =  I  T 

CONTINL"7 

GO  TO   426 

cTOP 

END 
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